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Preface 


The  work  described  in  this  report  was  authorized  as  part  of  the  Civil 
Works  Research  and  Development  Program  by  Headquarters,  U.S.  Army 
Corps  of  Engineers  (HQUSACE).  This  report  summarizes  wave  breaking 
investigations  performed  under  Work  Unit  31672,  “Nearshore  Waves  and 
Currents,*  at  the  Coastal  Engineering  Research  Center  (CERC)  of  foe  U.S. 
Army  Engineer  Waterways  Experiment  Station  (WES).  Messrs.  John  H. 
Lockhart,  Jr.,  John  G.  Housley,  Barry  W.  Holliday,  and  David  A.  Roellig 
were  HQUSACE  Technical  Monitors.  Ms.  Carolyn  M.  Holmes  was  foe 
CERC  Program  Manager. 
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Ms.  Jane  McKee  Smith,  Research  Hydraulic  Engineer,  CERC.  The  study 
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Charles  C.  Calhoun,  Jr.,  Director  and  Assistant  Director,  CERC,  respective¬ 
ly;  and  under  foe  direct  supervision  of  Mr.  H.  Lee  Butler,  Chief,  Research 
Division,  and  Mr.  Bruce  A.  Ebersole,  Chief,  Coastal  Processes  Branch, 
CERC.  Laboratory  measurements  of  waves  on  reefs  and  breaking  of  multiple 
wave  trains  were  made  by  Messrs.  John  M.  Heggins  and  Donald  L.  Ward, 
Wave  Research  Branch,  CERC,  respectively.  Ms.  Mary  A.  Cialone,  Mr.  W. 
Gray  Smith,  and  Mr.  Bruce  A.  Ebersole,  Coastal  Processes  Branch,  CERC, 
provided  technical  review  of  this  report. 

Dr.  Robert  W.  Whalin  was  Director  of  WES.  COL  Leonard  G.  Hassell, 
EN,  was  Commander. 


1  Introduction 


Background 

Waves  approaching  the  coast  increase  in  steepness  as  the  water  depth  de¬ 
creases.  When  the  wave  height  is  approximately  equal  to  the  water  depth,  the 
wave  breaks,  dissipating  wave  energy  and  inducing  longshore  currents  and  an 
increase  in  the  mean  water  level.  The  surf  zone  is  the  region  from  shoreline 
to  die  seaward  boundary  of  wave  breaking.  Within  the  surf  zone,  wave 
breaking  is  the  dominant  feature.  The  surf  zone  is  die  most  dynamic  coastal 
region  with  sediment  transport  and  bathymetry  change  driven  by  breaking 
waves  and  wave-induced  currents. 

Surf  zone  wave  conditions  are  required  for  estimating  potential  storm 
damage  (flooding  and  wave  damage),  calculating  shordine  evolution  and 
cross-shore  beach  profile  change,  designing  coastal  structures  (jetties,  groins, 
seawalls)  and  beach  fills,  and  developing  shordine  management  policy.  To 
provide  improved  estimates  of  surf  zone  waves,  the  Nearshore  Waves  and 
Currents  Work  Unit  in  the  Coastal  Flooding  and  Storm  Protection  Program 
developed,  validated,  and  improved  the  numerical  modeling  of  wave  breaking 
and  decay  through  laboratory  experiments  and  field  studies.  Areas  of  research 
included  incipient  breaker  indices,  surf  zone  wave  decay  expressions,  spectral 
shapes  in  the  surf  zone,  wave  breaking  on  reefs,  and  decay  of  multiple  wave 
trains. 


Scope 

The  purpose  of  this  report  is  to  summarize  the  nearshore  wave  breaking 
and  decay  research  performed  under  the  Nearshore  Waves  and  Currents  Work 
Unit.  Chapter  2  describes  incipient  wave  breaking  and  breaker  types.  Also, 
expressions  are  given  for  die  breaker  depth  index  and  breaker  height  index. 
Guidance  is  given  for  both  regular  and  irregular  waves.  The  decay  in  wave 
height  through  the  surf  zone  is  discussed  in  Chapter  3.  Regular  wave  methods 
include  the  simple  constant  ratio  of  wave  height  to  water  depth  and  die  energy 
flux  method.  Irregular  wave  methods  include  statistical  and  spectral  approach¬ 
es.  Chapter  4  applies  wave  height  decay  methodology  to  the  special  problem 
of  wave  breaking  on  reefs.  Reefs  have  steep  offshore  slopes  and  rough,  flat 
reef  tops,  so  the  wave  characteristics  differ  from  typical  sandy  beaches. 
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Chapter  S  presents  results  on  the  decay  of  multiple  wave  trains.  Generally, 
more  than  one  wave  train  exists  at  a  coastal  site,  but  no  guidance  exists  on 
how  to  apply  wave  breaking  to  multiple  wave  trains.  This  chapter  discusses 
laboratory  tests  on  the  decay  of  multiple  wave  trains. 
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2  Incipient  Wave  Breaking 


As  waves  approach  a  beach,  they  increase  in  steepness  (ratio  of  wave 
height  H  1  to  wavelength  L)  -  length  decreases  and  height  may  increase. 
Waves  break  as  they  reach  a  limiting  steepness  which  is  a  function  of  the 
relative  depth  (ratio  of  the  depth  d  to  the  wave  length)  and  the  beach  slope 
tan  /S.  Incipient  breaking  can  be  defined  several  ways  (Singamsetti  and  Wind 
1980): 

a.  The  point  where  the  wave  cannot  further  adapt  to  the  changing  bottom 
configuration  and  starts  to  disintegrate. 

b.  The  point  where  the  horizontal  component  of  the  water  particle  veloci¬ 
ty  at  the  crest  becomes  greater  than  the  wave  celerity. 

c.  The  point  where  the  wave  height  is  maximum. 

d.  The  point  where  part  of  the  wave  front  becomes  vertical. 

e.  The  point  where  the  radiation  stresses  start  to  decrease. 

/.  The  point  where  the  water  particle  acceleration  at  the  crest  tends  to 
separate  the  particles  from  the  water  surface. 

g.  The  point  where  the  pressure  at  the  free  surface  given  by  the  Bernoulli 
equation  is  incompatible  with  the  atmospheric  pressure. 

Determination  of  the  incipient  break  point  is  subjective  since  waves  generally 
break  gradually. 


Breaker  Type 

Breaker  type  refers  to  the  form  of  the  wave  at  breaking.  Wave  breaking 
may  be  classified  into  four  types  (Galvin  1968):  spilling,  plunging,  collaps¬ 
ing,  and  surging.  In  spilling  breakers  the  wave  crest  becomes  unstable  and 
cascades  down  the  shoreward  face  of  the  wave  producing  a  foamy  water  sur- 


1  For  convenience,  symbols  and  abbreviations  are  listed  in  the  notation  (Appendix  B). 
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ca'cades  down  the  shoreward  face  of  the  wave  producing  a  foamy  water  sur¬ 
face.  In  plunging  breakers  the  crest  curls  over  the  shoreward  face  of  the 
wave  and  falls  into  the  base  of  the  wave,  resulting  in  a  high  splash.  In  col¬ 
lapsing  breakers  the  crest  remains  unbroken  while  the  lower  part  of  the  shore¬ 
ward  face  steepens  and  then  foils,  producing  an  irregular  turbulent  water 
surface  In  surging  breakers  the  crest  remains  unbroken  and  the  front  face  of 
the  wave  advances  up  the  beach  with  minor  breaking. 

The  surf  similarity  parameter  £»,  defined  as 


where  the  subscript  o  denotes  deepwater  parameters,  is  used  to  distinguish 
breaker  type  (Galvin  1968,  Battjes  1974).  The  breaker  type  is  estimated  by 

surging/collapsing  for  >  3.3 

plunging  for  0.5  <  <  3.3  (2) 

spilling  for  <  0.5 

As  expressed  in  Equation  2,  spilling  breakers  tend  to  occur  for  high-steepness 
waves  on  low-slope  beaches.  Plunging  breakers  occur  on  steeper  beaches 
with  intermediately  steep  waves,  and  surging  breakers  occur  for  low  steepness 
waves  on  steep  beaches.  Extremely  low  steepness  waves  may  not  break,  but 
instead  reflect  from  the  beach,  forming  a  standing  wave. 

Spilling  breakers  differ  little  in  fluid  motion  from  unbroken  waves 
(Divoky,  LeMdhautg,  and  Lin  1970)  and  generate  less  turbulence  near  the 
bottom  and  thus  tend  to  be  less  effective  in  suspending  sediment  than  plunging 
or  collapsing  breakers.  The  most  intense  local  fluid  motions  are  produced  by 
a  plunging  breaker.  As  it  breaks,  the  crest  of  the  plunging  wave  acts  as  a 
free-falling  jet  that  scours  a  trough  into  the  bottom.  The  transition  from  one 
breaker  type  to  another  is  gradual  and  without  distinct  dividing  lines.  The 
direction  and  magnitude  of  the  local  wind  can  affect  breaker  type.  Onshore 
winds  cause  waves  to  break  in  deeper  water  depths  and  spill,  while  offshore 
winds  cause  waves  to  break  in  shallower  depths  and  plunge  (Galloway,  Col¬ 
lins,  and  Moran  1989;  Douglass  1990). 


Breaker  Index 


Many  studies  have  been  performed  to  develop  relationships  to  predict  the 
wave  height  at  incipient  breaking.  The  term  breaker  index  is  used  to  describe 
nondimensional  breaker  height.  Two  common  indices  are  of  the  form 


n  - 


<*> 


(3) 


Chapter  2  Incipient  Wave  Breaking 


where 


yb  -  breaker  depth  index 

Hk  =  wave  height  at  incipient  breaking 

dk  =  mean  water  depth  at  incipient  breaking 


0»- 


Hj, 


(4) 


in  which  is  the  breaker  height  index. 

Early  studies  on  breaker  indices  were  conducted  using  solitary  waves. 
McCowan  (1891)  theoretically  determined  the  breaker  depth  index  as 
yk  =  0.78  for  a  solitary  wave  traveling  over  a  horizontal  bottom.  Munk 
(1949)  derived  the  expression  Qb  =  0.3  (HJLJ,/3  for  the  breaker  height  index 
from  solitary  wave  theory.  Subsequent  studies,  based  on  periodic  waves,  by 
Iversen  (1952),  Goda  (1970),  Weggel  (1972),  Singamsetti  and  Wind  (1980), 
Sunamura  (1980),  Smith  and  Kraus  (1991),  and  others  have  established  that 
breaker  indices  depend  on  beach  slope  and  incident  wave  steepness. 

From  monochromatic  laboratory  data  on  smooth,  plane  slopes,  Weggel 
(1972)  derived  the  following  expression  for  the  maximum  breaker  depth  index 

Y  .  b  -  a  A.  (5) 

*  *  J* 


for  tan  /S  <.  1/10  and  HJL0  <,  0.06,  where  T  is  wave  period,  and  g  is  gravi¬ 
tational  acceleration.  The  parameters  a  and  b  were  empirically  determined 
functions  of  the  beach  slope,  given  by 

a  =  43.8  (1  -  e‘,9*“p) 


b  - 


1.56 

(1  ♦  e'IM  »  p) 


(7) 


The  breaking  wave  height  Hh  is  contained  on  both  sides  of  Equation  5,  so  the 
equation  must  be  solved  iteratively.  Figure  1  graphically  shows  the  breaker 
depth  index  dependency  on  wave  steepness  and  bottom  slope  (Equation  5). 

Smith  and  Kraus  (1991)  combined  monochromatic  laboratory  data  from 
Iversen  (1952),  Horikawa  and  Kuo  (1966),  Galvin  (1969),  Saeki  and  Sasaki 
(1973),  Iwagaki  et  al.  (1974),  Walker  (1974),  Singamsetti  and  Wind  (1980), 
Mizuguchi  (1980),  Maruyama  et  al.  (1983),  Visser  (1982),  and  Stive  (1985) 
with  data  collected  in  their  own  experiments  to  empirically  obtain 
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figure  1 .  Breaker  depth  index  as  a  function  of  HJ(g  T3)  (Equation  5) 
(Weggel  1972) 


Figure  2.  Breaker  depth  index  as  a  function  of  HJLC  (Equation  8)  (Smith  and 
Kraus  1991) 


Chapter  2  Incipient  Wave  Breaking 


for  0.0007  £  HJL0  £  0.0921  and  1/80  £  tan  /3  £  1/10.  Equation  8  is 
similar  in  form  to  Equation  S,  but  Hk  appears  only  in  yt  so  the  equation  can 
be  solved  without  iteration.  Equation  8  generally  gives  values  of  yt  that  are 
approximately  5  to  20  percent  smaller  than  Equation  5  for  1/80  £  tan  /9  £ 
1/10.  Figure  2  graphically  shows  the  breaker  depth  index  dependency  on 
wave  steepness  and  bottom  slope  (Equation  8). 

Komar  and  Gaughan  (1972)  derived  a  semi-empirical  relationship  for  the 
breaker  height  index  from  linear  wave  theory 

t 

5  (9) 


The  coefficient  0.56  was  determined  empirically  from  laboratory  and  field 
data.  Smith  and  Kraus  (1991)  empirically  determined  the  breaker  height  index 

fjj  \CM  -  o.«s  « 

Q„  =  (0.34  +  2.47  tan  0)  U-iJ  (10) 

based  on  the  same  monochromatic  laboratory  data  used  to  develop  Equation  8. 
Equations  9  and  10  give  similar  results  for  steep  beaches  and  low-steepness 
waves,  and  Equation  10  predicts  smaller  values  of  Q*  for  high-steepness  waves 
on  mild  slopes. 

In  irregular  seas,  incipient  breaking  may  occur  over  a  wide  zone  as  indi¬ 
vidual  waves  of  different  heights  and  periods  reach  their  steepness  limits.  In 
the  saturated  breaking  zone  for  irregular  waves  (the  zone  where  essentially  all 
waves  are  breaking),  the  wave  height  may  be  related  to  the  local  depth 

*  0.42  d  (11) 


for  root-mean-square  wave  height  (Thornton  and  Guza  1983)  or,  approximate¬ 
ly. 

=  0.6  d  (12) 


for  zero-moment  wave  height.  Davis,  Smith,  and  Vincent  (1991)  determined 
that  wave  breaking  occurred  for  height  to  sttll-water  depth  ratios  exceeding 


0.243 


\-o.u 


8  T, 


/  / 


(13) 
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where  hk  is  the  still-water  depth  at  breaking,  h  is  the  still-water  depth,  and  Tf 
is  the  peak  wave  period,  based  on  irregular  wave  laboratory  data  with  tan  /S  = 
1/30  and  0.008  £  Qi,JJLe  £  0.044.  Some  variability  in  and 
with  wave  steepness  and  beach  slope  should  be  expected;  however,  no  com¬ 
prehensive  study  has  been  performed. 
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3  Surf  Zone  Wave 
Transformation 


Following  incipient  wave  breaking,  the  wave  shape  changes  rapidly  to 
resemble  a  bore  (Svendsen  1984).  The  wave  profile  becomes  sawtooth  in 
shape  with  the  leading  edge  of  die  wave  crest  becoming  nearly  vertical  (Fig¬ 
ure  3).  The  wave  may  continue  to  dissipate  energy  to  die  shoreline  or,  if  the 
water  depth  again  increases  as  in  die  case  of  a  barred  beach  profile,  the  wave 
may  cease  breaking,  re-form,  and  break  again  on  the  shore.  The  transforma¬ 
tion  of  the  wave  height  through  the  surf  zone  impacts  wave  setup,  runup, 
nearshore  currents,  and  sediment  transport. 


Regular  Waves 

Similarity  method 

The  simplest  method  for  predicting  wave  height  through  the  surf  zone,  an 
extension  of  Equation  3  shoreward  of  incipient  breaking  conditions,  is  to 
assume  a  constant  height-to-depth  ratio  from  the  break  point  to  the  shore 

H=ybd  (14) 

This  method  has  been  used  successfully  by  Longuet-Higgins  and  Stewart 
(1963)  to  calculate  setup,  and  by  Bowen  (1969),  Longuet-Higgins  (1970a,b), 
and  Thornton  (1970)  to  calculate  longshore  currents.  The  similarity  method  is 
applicable  only  for  monotonically  decreasing  water  depth  through  the  surf 
zone  and  gives  best  results  for  a  beach  slope  of  approximately  1/30.  On 
steeper  slopes.  Equation  14  tends  to  underestimate  the  wave  height,  and  on 
shallower  slopes  or  barred  topography,  it  tends  to  overestimate  the  wave 
height.  Equation  14  is  based  on  the  assumption  that  the  wave  height  is  zero  at 
the  mean  shoreline.  Camfield  (1991)  shows  that  a  conservative  estimate  of 
die  wave  height  at  the  still-water  shoreline  is  0.20  Hh  for  1/100  £  tan  0  £ 
1/10. 
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Smith  and  Kraus  (1988)  developed  an  improved  analytical  expression  for 
wave  height  in  the  surf  zone  for  planar  beaches  that  includes  the  effect  of 
beach  slope,  which  is  neglected  in  Equation  14 


H  - 


(IS) 


where 


n 


0.66  yb  + 


0-043  y  b 
tan  0 


+  0.032 

tan  $ 


(16) 


Equations  IS  and  16  require  that  the  incipient  breaking  conditions  (Hb,  hh,  and 
yd  be  specified  (e.g..  Equation  8). 


Energy  flux  method 

A  more  general  method  for  predicting  the  wave  height  through  die  surf 
zone  is  solution  of  the  steady-state  energy  balance  equation 

d(E  C  cos  0)  , 

*  e  — O  V 

dx 


where 

E  =  wave  energy  per  unit  surface  area 
Cf  =  wave  group  speed 

B  =  wave  direction  relative  to  shore-normal 

jc  =  cross-shore  coordinate 

6  =  energy  dissipation  rate  per  unit  surface  area  due  to  wave  breaking 

The  wave  energy  flux  E  Ct  may  be  specified  from  linear  or  higher  order  wave 
theory.  LeMdhautd  (1962)  approximated  a  breaking  wave  as  a  hydraulic  jump 
and  substituted  the  dissipation  of  a  hydraulic  jump  for  6  in  Equation  17  (see 
also  Divoky,  LeMdhautd,  and  Lin  1970;  Hwang  and  Divoky  1970;  Svendsen, 
Madsen,  and  Hansen  1978). 

Dally,  Dean,  and  Dalrymple  (1985)  define  the  dissipation  rate  as 

6  -  -  (EC.  -  EC,  J  (18) 

d  *  *’ 


where  x  is  an  empirical  decay  coefficient,  found  to  have  the  value  0.15,  and 
E  CtJ  is  the  stable  energy  flux  associated  with  the  stable  wave  height 


(19) 
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where  T  is  an  empirical  coefficient  with  a  value  of  approximately  0.40.  The 
stable  wave  height  is  the  height  at  which  a  wave  stops  breaking  and  re-forms. 
This  approach  is  based  on  the  assumption  that  energy  dissipation  is  propor¬ 
tional  to  the  difference  between  the  local  energy  flux  and  the  stable  energy 
flux.  Applying  linear,  shallow-water  theory  and  assuming  normal  wave  inci¬ 
dence,  die  Dally,  Dean,  and  Dalrymple  model  is  given  by 

(H*  -  r=  d*)  for  (20) 

=  0  for  H  < 

This  approach  has  been  successful  in  modeling  wave  transformation  over 
irregular  beach  profiles,  including  bars  (e.g.,  Ebersole  1987,  Larson  and 
Kraus  1991,  Dally  1992). 


Irregular  Waves 

Transformation  of  irregular  waves  through  the  surf  zone  may  be  analyzed 
or  modeled  with  either  a  statistical  (individual  wave  or  wave  height  distribu¬ 
tion)  or  a  spectral  (parametric  spectral  shape)  approach. 


Statistical  approaches 

Individual  waves.  The  most  straightforward  statistical  approach  is  the 
transformation  of  individual  waves  through  the  surf  zone.  Individual  waves 
seaward  of  breaking  may  be  measured  direcdy,  randomly  chosen  from  a 
Rayleigh  distribution,  or  chosen  to  represent  wave  height  classes  in  the  Ray¬ 
leigh  distribution.  Then  the  individual  waves  are  independently  transformed 
through  the  surf  zone  using  Equation  17.  The  distribution  of  wave  heights 
can  be  calculated  at  any  point  across  the  surf  zone  by  recombining  the  individ¬ 
ual  waves  into  a  distribution  to  calculate  wave  height  statistics  (e.g.,  average 
of  the  highest  1/10  waves  H,nth  significant  wave  Hw,  and  This  method 
does  not  make  any  a  priori  assumptions  about  the  wave  height  distribution  in 
the  surf  zone.  The  individual  wave  method  has  been  applied  and  validated 
with  field  data  by  Dally  (1990,  1992),  and  Larson  and  Kraus  (1991). 

A  numerical  model  called  NMLONG  (Numerical  Model  of  the  LONGshore 
current)  (Larson  and  Kraus  1991,  Kraus  and  Larson  1991)  is  available  for 
calculating  wave  breaking  and  decay  by  the  individual  wave  approach  apply¬ 
ing  the  Dally,  Dean,  and  Dalrymple  (1985)  wave  decay  model  (monochromat¬ 
ic  or  irregular  waves).  The  main  assumption  underlying  the  model  is  unifor¬ 
mity  of  waves  and  bathymetry  alongshore,  but  the  beach  profile  can  be  irregu¬ 
lar  across  the  shore  (e.g.,  longshore  bars  and  nonuniform  slopes).  The  model 
runs  on  a  personal  computer  and  has  a  convenient  menu  and  graphical  inter¬ 
face.  NMLONG  calculates  both  wave  transformation  and  longshore  current 
for  arbitrary  offshore  (input)  wave  conditions,  and  it  provides  a  plot  of  the 
results.  Figures  4  and  5  give  example  NMLONG  calculations  and  compari- 


12 


Chapter  3  Surf  Zona  Wave  Transformation 


0  20  40  60  80  100  120  140 

Distance  Offshore  (m) 


Figure  4.  Example  NMLONG  simulation  of  wave  height  transfor¬ 
mation,  Leadbetter  Beach,  California,  3  Feb  1980 
(Thornton  and  Guza  1 986)) 


Figure  5.  Example  NMLONG  simulation  of  wave  height  transform 
ation,  Duck,  North  Carolina,  14  Oct  1990 
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sons  to  field  measurements  of  wave  breaking  at  Leadbetter  Beach,  California 
(Thornton  and  Guza  1986)  and  Duck,  North  Carolina,  respectively.  The 
water  depth  at  Leadbetter  Beach  decreases  monotonically  to  the  beach.  At 
Duck,  the  beach  has  a  bar  and  trough  profile,  so  Figure  5  shows  two  zones  of 
wave  breaking  with  wave  re-formation  between  them.  As  with  all  numerical 
models,  proper  use  of  NMLONG  requires  careful  examination  of  the  related 
documentation. 

Wave  height  distribution.  A  second  statistical  approach  is  based  on  the 
assumption  of  a  known  wave  height  distribution  in  the  surf  zone.  The  Ray¬ 
leigh  distribution  is  a  reliable  measure  of  the  wave  height  distribution  in  deep 
water  and  at  finite  depths.  But,  as  waves  approach  breaking,  the  distribution 
of  wave  heights  shows  a  definite  weighting  toward  the  higher  end  of  the  distri¬ 
bution.  Also,  depth-induced  breaking  acts  to  limit  the  highest  wave  in  the 
distribution,  contrary  to  the  Rayleigh  distribution,  which  is  unbounded.  The 
surf  zone  wave  height  distribution  has  generally  been  represented  as  a  truncat¬ 
ed  Rayleigh  distribution  (e.g.,  Collins  1970,  Battjes  1972,  Kuo  and  Kuo  1974, 
Goda  1975).  Battjes  and  Janssen  (1978)  and  Thornton  and  Guza  (1983)  base 
the  distribution  of  wave  heights  at  any  point  in  the  surf  zone  on  a  Rayleigh 
distribution  or  a  truncated  Rayleigh  distribution  (truncated  above  a  maximum 
wave  height  for  the  given  water  depth).  A  percentage  of  the  waves  in  the 
distribution  is  designated  as  broken  (based  on  the  ratio  Hid),  and  energy 
dissipation  from  these  broken  waves  is  calculated  from  Equation  17  through  a 
model  of  dissipation  similar  to  a  periodic  bore.  The  dissipation  across  the 
surf  zone  specifies  H^,,  which  in  turn  is  used  to  parameterize  the  Rayleigh 
distribution.  This  method  has  been  validated  with  laboratory  and  field  data 
(Battjes  and  Janssen  1978,  Thornton  and  Guza  1983). 


Spectral  approaches 

In  shallow  water,  the  shape  of  the  wave  spectrum  is  influenced  by  nonlin¬ 
ear  transfers  of  wave  energy  from  the  peak  frequency  to  higher  and  lower 
frequencies  (Freilich  and  Guza  1984;  Freilich,  Guza,  and  Elgar  1990).  High¬ 
er  harmonics  (energy  peaks  at  multiples  of  the  peak  frequency)  appear  in  the 
spectrum  as  well  as  a  general  increase  in  the  energy  level  above  the  peak 
frequency  as  illustrated  in  Figure  6.  Figure  6  shows  wave  spectra  measured 
at  three  locations  in  a  large  wave  flume  with  a  sloping  sand  beach.  The  solid 
curve  is  the  incident  spectrum  ( d  =  3.0  m),  the  dotted  curve  is  the  spectrum 
at  the  zone  of  incipient  breaking  (d  =  1.7  m),  and  the  dashed  curve  is  within 
the  surf  zone  (d  =  1.4  m).  Measured  surf-zone  spectra  show  that  the  energy 
dissipated  at  each  frequency  is  roughly  proportional  to  the  energy  content  at 
that  frequency.  Presently,  no  formulation  is  available  for  the  dissipation  rate 
based  or.  spectral  parameters  for  use  in  Equation  17.  Therefore,  the  energy  in 
the  spectrum  is  often  limited  using  a  depth-limited  wave  height  (e.g., 

Equation  12). 

Davis,  Smith,  and  Vincent  (1991)  observed  that  laboratory  wave  spectra 
transformed  to  a  uniform  spectral  shape  within  the  surf  zone  for  a  given  peak 
spectral  period,  irrespective  of  the  total  spectral  energy  in  deep  water.  They 
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Figure  6.  Shallow-water  wave  spectra  (solid  -incident,  d  «  3.0  m;  dot  - 
incipient  breaking,  d  *  1.7  m;  dash  -  surf  zone,  d  *  1 .4  m) 

parameterized  the  spectrum  in  the  surf  zone  based  on  a  modification  of  the 
FRF  spectrum  (named  for  the  Field  Research  Facility  (FRF)  of  the  U.S. 
Army  Engineer  Waterways  Experiment  Station  Coastal  Engineering  Research 
Center)  (Miller  and  Vincent  1990).  The  FRF  spectrum  is  expressed  as 


UFXF 


(2D 


where 


EfXF  ~  spectral  energy  at  frequency/ 
a,  =  equilibrium  range  constant  (=  0.0029) 
spectral  wave  energy  parameter 
k  =  wave  number 

a  -  radian  wave  frequency 

/  =  frequency 

fp  =  peak  frequency 

7  =  spectral  peakedness  factor 

o  =  spectral  width  parameter 

Regression  analysis  was  used  to  correlate  the  spectral  parameters  0 W»  7-  °) 
to  the  incident  wave  characteristics.  In  application  of  the  FRF  spectrum  to 
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prototype  conditions,  the  parameter  is  equivalent  to  the  10-m  elevation 
wind  speed.  In  the  laboratory  experiments,  however,  the  wave  spectra  were 
generated  by  a  wave  paddle  and  not  the  wind.  Under  such  conditions,  /W  is 
a  spectral  parameter  representing  the  wave  energy,  but  not  related  to  wind 
speed.  The  regression  analysis  gave  the  following  results  for  the  limiting 
(maximum)  energy  levels  for  the  FRF  spectrum: 

¥  ■  -  W 

=  32 

where  C  is  wave  celerity  (=  L/T^ 


outside  surf  zone  (22) 

inside  surf  zone 


and 


o  e  0.11 


(24) 


for  0.008  ^  (H'JJL,  ^  0.044  and  tan  0  =  1/30.  Equations  23  and  24  have 
not  been  tested  for  a  wide  range  of  incident  spectral  shapes;  thus,  they  should 
be  used  with  caution.  Equations  21  through  24  are  applied  in  the  spectral 
wave  propagation  model  STWAVE  (Cialone  et  al.  1992)  to  limit  wave  energy 
due  to  breaking. 
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4  Wave  Attenuation 
Over  Reefs 


Many  tropical  coastal  regions  are  fronted  by  coral  reefs.  These  reefs  offer 
protection  to  the  coast  because  waves  break  on  die  reefs,  so  the  waves  reach¬ 
ing  the  shore  are  much  reduced  in  height.  Reefs  typically  have  steep  seaward 
slopes  with  broad,  flat  reef  tops  and  a  deeper  lagoon  shoreward  of  die  reef. 
The  transformation  of  waves  across  coral  reefs  is  a  complex  problem,  includ¬ 
ing  the  processes  of  refraction,  shoaling,  breaking,  energy  dissipation  by 
bottom  friction,  and  reflection. 


Wave  Processes 

As  waves  pass  from  deep  water  over  a  steep  reef  face  onto  die  reef  flat, 
the  waves  become  highly  nonlinear.  Wave  energy  is  dissipated  due  to  break¬ 
ing,  but  energy  is  also  transferred  to  both  higher  and  lower  frequencies  in  the 
wave  spectrum,  and  the  spectral  shape  becomes  flat  (Young  1989,  Hardy  and 
Young  1991).  The  peak  wave  period  shoreward  of  the  reef  face  may  become 
shorter  as  higher  harmonics  are  transmitted  as  free  waves  (Lee  and  Black 
1978),  or  the  period  may  increase  as  surf  beat  dominates  the  spectrum. 
Brealdng  waves  induce  a  setup  of  the  water  surface  over  the  reef,  and  differ¬ 
ences  in  breaking  characteristics  along  the  reef  can  cause  variations  in  wave 
setup,  producing  significant  longshore  currents.  Although  it  may  seem  that 
wave  reflection  off  a  nearly  vertical  reef  would  be  significant,  field  data 
(Young  1989;  Roberts,  Murray,  and  Suhayda  1975)  have  shown  the  reflected 
wave  height  to  be  on  the  order  of  only  10  percent  of  the  incident  height  (due 
to  the  porosity  of  the  reef).  Energy  losses  due  to  bottom  friction  are  usually 
negligible  in  wave  transformation  across  sandy  beach  profiles,  but  may  be 
significant  over  shallow,  rough  reef  flats  where  the  bottom  friction  coefficient 
may  be  an  order  of  magnitude  larger  than  for  a  sandy  bed  (Roberts,  Murray, 
and  Suhayda  1975;  Gerritsen  1980). 


Wave  breaking 

For  engineering  purposes,  the  most  significant  wave  transformation  process 
on  a  reef  is  generally  depth-limited  breaking.  Design  wave  heights  for  a 
breaking  wave  on  a  structure  are  often  determined  from  a  bottom  slope-depen- 
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dent  maximum  height-to-depth  ratio  at  the  toe  of  the  structure  (or  just  o' 
shore).  Over  a  flat  reef,  this  would  predict  breaking  wave  heights  of  0.78 
times  the  depth.  This  breaking  wave  height  ratio  is  overly  conservative  for 
design  wave  heights  on  the  shoreward  edge  of  a  wide  reef  or  in  the  lagoon 
behind. 

The  concept  of  a  constant  height-to-depth  ratio  in  the  surf  zone  is  in¬ 
correct;  prototype  data  show  the  height-to-depth  ratio  varying  between  1.1  and 
0.4  across  reeft  (Gerritsen  1980,  Hardy  et  al.  1990).  Similar  to  die  case  of 
waves  breaking  on  a  barred  beach,  waves  on  a  reef  flat  will  break,  dissipating 
energy,  and  then  re-form  as  they  travel  across  die  reef.  Wave  height  will 
decay  quickly  on  the  outer  portion  of  the  reef  until  it  reaches  a  stable  value. 
On  die  inner  portion  of  the  reef,  the  re-formed  wave  will  decay  slowly  due  to 
bottom  friction.  The  breaking  and  re-formation  process  is  strongly  dependent 
on  the  width  of  the  reef  and  the  water  depth  over  the  reef.  To  accurately 
estimate  wave  heights  on  the  reef  or  in  a  lagoon,  it  is  necessary  to  model 
transformation  across  the  entire  reef  and  to  represent  wave  setup  (driven  by 
the  gradient  in  wave  height).  Wave  height  estimates  based  only  on  incident 
wave  conditions  and  still-water  depths  over  the  reef  will  not  be  reliable  across 
the  entire  reef. 


Breaker  models 

Wave  breaking  and  re-formation  on  a  reef  are  similar  to  the  process  on  a 
barred  beach  (see  Chapter  3).  Gerritsen  (1980)  first  applied  wave  breaking 
methods  developed  for  mildly  sloping  beaches  to  reefs,  with  the  inclusion  of 
dissipation  due  to  bottom  friction.  Gerritsen  applied  the  random  breaking 
wave  model  developed  by  Battjes  and  Janssen  (1978),  but  found  that  the  trun¬ 
cated  Rayleigh  distribution  of  wave  heights  assumed  by  Battjes  and  Janssen 
was  a  poor  representation  of  broken  waves  over  a  reef.  Young  (1989)  used  a 
similar  approach,  but  included  a  check  to  "turn  off  wave  breaking  when  the 
height-to-depth  ratio  was  less  than  0.78,  to  simulate  wave  re-formation. 

The  Dally,  Dean,  and  Dalrymple  (1985)  wave  breaking  and  re-formation 
model  has  been  extensively  verified  for  plane  beaches,  composite  beach 
slopes,  and  barred  beach  profiles.  The  model  is  based  on  the  transformation 
of  individual  waves,  but  may  be  applied  to  a  random  wave  field  using  a  statis¬ 
tical  approach.  The  statistical  approach  requires  specification  of  the  wave 
height  distribution  in  the  offshore  region,  but  does  not  impose  a  specified 
distribution  in  the  surf  zone.  The  advantages  of  die  Dally,  Dean,  and 
Dalrymple  model  are:  a)  extensive  verification  for  a  variety  of  beach  config¬ 
urations,  b)  no  a  priori  specification  of  the  wave  height  distribution  in  the 
surf  zone,  and  c)  the  individual  wave  approach  allows  calculation  of  the  wave 
height  distribution  and  statistical  wave  height  parameters  (/fw,  Hm,  in 
the  surf  zone.  Due  to  these  advantages,  the  Dally,  Dean,  and  Dalrymple 
model  is  recommended  to  calculate  wave  attenuation  over  reefs. 
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Modeling  Methodology 


The  methodology  described  here  includes  the  processes  of  wave  shoaling, 
refraction,  depth-limited  breaking,  and  bottom  friction.  The  assumptions 
include:  linear  wave  theory,  steady-state  wave  conditions,  Rayleigh  wave 
height  distribution  in  the  offshore,  and  longshore  homogeneity.  The  method 
neglects  energy  shifts  within  the  wave  spectrum,  wave-current  interaction,  and 
wave  reflection  and  scattering. 


The  steady-state  energy  balance  equation  governing  wave  propagation  is 
given  by 


d(E  Cg  cos  6) 
dx 


P  Cf  /  2  it  H  V 

6  it  (Tsinh  (k  d)) 


(25) 


where 

E  —  wave  energy  (=  p  g  If/8  for  linear  wave  theory) 

Cg,  =  wave  group  speed  associated  with  the  stable  wave  height 

p  =  density  of  water 

Cf  =  bottom  friction  coefficient 

The  first  term  on  the  right  side  of  Equation  25  is  the  energy  dissipation  due  to 
wave  breaking  (Dally,  Dean,  and  Dalrymple  1985)  (see  Equations  17  and  18), 
and  the  second  term  is  energy  dissipation  due  to  bottom  friction  (Gerritsen 
1980,  Thornton  and  Guza  1983).  The  application  of  Equation  25  to  calculate 
random  wave  transformation  across  reefs  is  based  on  the  approach  of  Larson 
and  Kraus  (1991).  The  input  parameters  required  include  the  cross-shore 
profile  of  the  reef  and  the  offshore  wave  period,  mean  direction,  and  root- 
mean-square  height  H",.  The  wave-breaking  parameters  (height-to-depth  ratio 
for  incipient  breaking  y*,  *,  and  T)  and  bottom  friction  coefficient  must  also 
be  specified  (Gerritsen  suggests  values  of  Q  =  0.05  to  0.25  for  coral  reefs). 
From  the  specified  offshore  H^,  a  Rayleigh  distribution  of  wave  heights  is 
determined.  Individual  wave  heights  are  randomly  chosen  from  the  Rayleigh 
distribution,  and  each  individual  wave  is  transformed  independently.  Wave 
angles  are  calculated  by  Snell’s  law,  wave  setup  is  calculated  from  the  cross- 
shore  balance  of  momentum,  driven  by  cross-shore  gradients  in  wave  height 
(Longuet-Higgins  and  Stewart  1964),  and  wave  height  is  calculated  from 
Equation  25.  The  wave  height  statistic  is  determined  across  the  reef  by 
combining  the  transformed  individual  waves.  Other  wave  height  statistics, 
e.g.,  Hm  or  Hmo,  may  also  be  calculated.  Generally  100  or  more  individual 
waves  are  required  for  stable  mean  statistics. 


Sample  Results 

Limited  validation  of  the  method  described  above  was  performed  using 
laboratory  data  from  a  flume  study  with  a  configuration  replicating  the  reef  at 
Agat,  Guam,  and  field  data  from  Yonge  Reef,  Australia  (Young  1989).  The 
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laboratory  study  was  conducted  in  an  46-cm-wide  flume  consisting  of  21.3  m 
of  flat  bottom,  1.5  m  of  1/5  slope,  10.2  m  of  1/30  slope,  12.9  m  of  reef  flat, 
and  3.3  m  (covered  with  wave  absorber)  of  1/30  slope.  The  water  depth  in 
the  deepest  portion  of  the  flume  was  0.69  to  0.64  m  with  a  depth  on  the  reef 
flat  of  0.05  to  0.005  m.  Wave  periods  ranged  from  1.1  to  2.5  sec  and  heights 
ranged  from  0.02  to  0. 12  m.  Laboratory  data  are  summarized  in 
Appendix  A.  The  model  was  applied  with  the  standard  breaking  parameters 
(*  =  0. 15  and  T  =  0.4)  and  an  incipient  breaker  index  y*  =  1.0  (a  breaking 
index  of  0.8  is  commonly  used  on  gently  sloping  beaches,  but  an  index  of  1.0 
is  more  appropriate  for  steeper  seaward  reef  slopes).  The  bottom  friction 
coefficient  was  set  to  0.01,  which  is  a  typical  value  for  smooth  slopes  such  as 
the  lab  configuration.  Figures  7  through  9  show  selected  results  for  the  labo¬ 
ratory  data.  The  agreement  between  laboratory  measurements  and  model 
results  is  excellent.  The  solid  line  is  the  modeled  wave  height,  the  symbols 
are  the  measured  wave  height,  the  chain-dot  line  is  the  modeled  setup,  and  the 
dotted  line  is  the  still-water  level.  These  results  are  typical  for  water  depths 
greater  than  3  cm  on  the  reef  flat.  For  shallow  water  depths,  the  model  un¬ 
derpredicted  the  (small)  measured  wave  heights.  For  very  shallow  depths, 
wave  energy  at  the  incident  frequency  is  almost  entirely  dissipated  and  low- 
frequency  energy  (which  is  not  included  in  the  model)  dominates.  Figure  10 
is  a  scatter  plot  of  the  calculated  versus  measured  results  for  69  laboratory 
tests  (345  data  points).  The  figure  shows  the  good  correlation  for  the  higher 
wave  heights  (depths  greater  than  3  cm)  and  the  underprediction  of  the  low 
wave  heights  (depths  less  than  3  cm). 

Figure  1 1  shows  results  from  a  field  experiment  conducted  on  Yonge  Reef, 
which  is  part  of  the  Great  Barrier  Reef  in  Australia.  This  is  one  of  four  cases 
reported  by  Young  (1989).  The  incident  H w  =  2.05  m  and  the  depth  over 
the  reef  was  1.05  m.  The  wave  measurement  was  taken  in  the  lagoon  on  the 
leeward  side  of  the  reef.  The  wave  breaking  parameters  used  in  the  model 
were  identical  to  those  applied  for  the  laboratory  cases,  and  the  bottom  fric¬ 
tion  coefficient  was  0.05,  which  is  equivalent  to  the  value  suggested  by 
Young.  As  in  the  case  of  the  laboratory  results,  the  model  compares  well 
with  the  measurements. 


Recommendations 


For  engineering  purposes,  the  breaking  wave  model  of  Dally,  Dean,  and 
Dalrymple  can  be  used  to  calculate  the  attenuation  of  waves  over  reefs.  For 
very  small  water  depths  over  the  reef,  the  model  may  underpredict  wave 
height  as  nonlinear  processes  dominate.  The  inclusion  of  bottom  friction  in 
the  energy  balance  equation  improves  estimates  of  wave  height  across  the  reef 
flat,  but  may  not  be  critical  for  engineering  application.  There  are  insufficient 
measurements  to  give  general  guidance  of  bottom  friction  coefficients,  so  site- 
specific  field  measurements  are  recommended  to  determine  Cp 
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Figure  9.  Laboratory  test  37,  Agat,  Guam 


Figure  10.  Calculated  versus  measured  laboratory  results 
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Figure  1 1 .  Field  test.  August  1 985,  Yonge  Reef,  Australia  fYoung  1 989) 
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5  Shoaling  and  Decay 
of  Multiple  Wave  Trains 


Little  information  is  available  on  the  refraction,  shoaling,  and  breaking  of 
multiple-wave-train  systems.  Observations  in  relatively  shallow  water  (4-8  m) 
indicate  that  approximately  65  percent  of  wave  fields  are  comprised  of  two  or 
more  distinct  wave  trains  well  separated  in  the  frequency  domain  (Thompson 
1980).  Thompson’s  statistics  were  from  representative  sites  on  the  Atlantic, 
Pacific,  and  Gulf  of  Mexico  coasts  of  the  United  States,  and  may  be  indicative 
of  open  ocean  and  coastal  conditions  worldwide. 

Several  approaches  have  been  used  to  estimate  the  transformation  of  multi¬ 
ple-peaked  wave  systems.  The  traditional,  pervasive  approach  is  to  assign  all 
wave  energy  to  the  dominant  peak,  then  refract  and  shoal  this  representative 
wave,  applying  a  breaking  criterion  where  needed  (Larson,  Kraus,  and  Byrnes 
1989).  Another  method  is  to  divide  the  wave  field  into  two  (or  more)  wave 
trains  and  treat  each  one  separately  (independent  and  noninteracting),  with  a 
breaking  criterion  applied  to  each  wave  train  independently  and  the  results 
superimposed  (Gravens,  Scheffner,  and  Hubertz  1989)  or  to  consider  only  the 
most  significant  peak  and  disregard  other  peaks  (Kraus  et  al.  1988).  Some 
existing  spectral  wave  models  compute  each  frequency  and  direction  band  of 
the  spectrum  independently  (O’Reilly  and  Guza  1991,  Resio  1988).  Breaking 
may  be  treated  with  a  limiting  spectral  form,  an  energy  sink  term  in  the 
energy  balance  equation,  or  a  simple  renormalization  of  the  spectrum  to  match 
an  equivalent  monochromatic  wave  of  depth-limited  height. 


Laboratory  Experiments 

Laboratory  tests  were  conducted  to  examine  the  shoaling  and  decay  of 
multiple  wave  trains  in  a  wave  flume  (Smith  and  Vincent  1992)  and  a  wave 
basin  (Briggs,  Smith,  and  Green  1991).  Double-peaked  spectra  were  deter¬ 
mined  by  superimposing  two  spectra  of  the  TMA  form  (Bouws  et  al.  1985). 
Two  wave  period  combinations  were  used,  Tf  =  2.5/1.25  sec  and  Tf  — 
2.5/1.75  sec,  with  two  total  zero-moment  wave  heights,  H „  =  15.2  cm  and 
9.2  cm.  The  relative  splits  of  energy  between  the  low-frequency  and  high- 
frequency  peaks  were  l/3:2/3,  1/2: 1/2,  and  2/3: 1/3.  In  die  basin  tests, 
directional  spreading  of  0  deg,  20  deg,  and  40  deg  and  mean  wave  directions 
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of  0  deg  and  20  deg  were  included  in  the  test  series.  Both  the  flume  and 
basin  experiments  had  similar  bottom  configurations,  with  a  flat  concrete 
bottom  near  the  generators  and  a  1/30  concrete  slope.  Wave  gauges  were 
placed  across  the  slopes  to  measure  the  wave  transformation  with  water  depth. 
In  the  wave  basin,  two  longshore  arrays  (one  on  the  flat  and  one  on  the  slope) 
were  also  deployed  to  measure  wave  direction. 


Wave  Spectra 

Figures  12  through  15  summarize  the  spectral  transformation  for  four  of 
the  cases  simulated  in  the  laboratory  flume.  Figures  12  and  13  present  cases 
corresponding  to  dual  wave  trains  with  peak  periods  of  2.5  sec  and  1.25  sec 
with  —  15.2  cm.  Similarly,  Figures  14  and  15  present  cases  with  dual 
wave  trains  with  peak  periods  of  2.5  sec  and  1.75  sec  and  H^,  =  9.2  cm.  In 
Figures  13  and  15  approximately  two-thirds  of  the  energy  is  initially  in  the 
low  frequency  peak,  and  in  Figures  12  and  14  approximately  one-third  of  the 
energy  is  initially  in  the  low  frequency  peak.  In  each  figure,  only  four  of  the 
measured  spectra  are  overlaid  —  the  deepest  measurement  ( h  —  61  cm),  die 
shallowest  measurement  (A  =  6.1  cm),  and  two  intermediate-depth  measure¬ 
ments  (h  —  18.3  cm  and  9.1  cm)  —  but  these  are  sufficient  to  illustrate  the 
basic  trends  seen  in  the  complete  set  of  data. 


Figure  12.  Energy  spectra  for  T0  *=  2.5/1.25  sec,  H„,  *  15.2  cm  (/»  « 

61  cm  (solid),  18.3  cm  (dash),  9.1  cm  (dot),  6.1  cm  (chain-dot)) 
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Figure  13.  Energy  spectra  for  7#  *  2.5/1.25  sec.  «  15.2  cm  (h  m 

61  cm  (solid),18.3  cm  (dash),  9.1  cm  (dot),  6.1  cm  (chain-dot)) 


Figure  14.  Energy  spectra  for  7#  -  2.5/1.75  sec,  H «  9.2  cm  (/»  - 

61  cm  (solid),  18.3  cm  (dash),  9.1  cm  (dot).  6.1  cm  (chain-dot)) 
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Figure  1 5.  Energy  spectra  for  T0  =  2.5/1 .75  sec,  =  9.2  cm  (h  «= 

61  cm  (solid),  18.3  cm  (dash),  9.1  cm  (dot),  6.1  cm  (chain-dot)) 

The  trend  in  the  spectral  transformation  is  very  similar  for  each  case. 
Between  the  depths  of  61.0  cm  and  18.3  cm,  there  is  a  general  increase  in 
energy  at  frequencies  both  above  and  below  the  primary  peaks.  Also,  distinct 
energy  peaks  appear  at  harmonic  frequencies  due  to  nonlinear,  near-resonance 
triad  interaction  (Freilich  and  Guza  1984).  For  example,  in  Figure  12  a 
secondary  peak  appears  at  1.2  Hz,  which  is  the  sum  frequency  of  the  two 
primary  peaks  at  0.4  Hz  and  0.8  Hz.  In  Figure  IS,  secondary  peaks  appear  at 
0.8  Hz  (self-interaction  of  0.4-Hz  peak),  0.97  Hz  (interaction  of  0.4-Hz  and 
0.57-Hz  peaks),  and  1.14  Hz  (self-interaction  of  0.4-Hz  peak).  These  ex¬ 
amples  do  not  show  energy  transfer  into  the  valley  between  the  primary  peaks 
as  shown  by  Elgar  and  Guza  (1985),  probably  because  harmonics  of  the  pri¬ 
mary  peaks  do  not  exist  in  this  valley.  For  Figures  12  and  13,  the  second 
peak  (at  0.8  Hz)  was  at  a  harmonic  of  the  first  peak  (at  0.4  Hz),  but  this  did 
not  seem  to  affect  the  spectral  transformation  since  Figures  14  and  15  show 
the  same  trends. 

Between  depths  of  18.3  cm  and  6.1  cm,  wave  breaking  occurred.  Wave 
breaking  lowered  the  general  energy  level  from  the  low-frequency  peak 
through  all  the  higher  frequencies  and,  by  the  shallowest  gauge,  eliminated  all 
energy  peaks  above  the  low-frequency  peak.  The  low-frequency  peak  became 
die  dominant  peak  of  the  spectrum  as  the  waves  shoaled  and  broke.  In  all  test 
cases,  only  the  low-frequency  peak  is  clearly  dominant  in  the  shallowest  meas¬ 
urement.  This  trend  was  evident  whether  the  peaks  were  closely  placed,  the 
low-frequency  peak  was  initially  the  lower-energy  peak  in  the  spectrum,  the 
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total  significant  wave  height  was  high  or  low,  directional  spreading  was  ap¬ 
plied,  or  the  peaks  had  the  same  or  different  mean  directions. 

Observations  show  that  the  dominance  of  the  low-frequency  peak  is  due 
primarily  to  the  loss  of  energy  in  the  high-frequency  peak.  Simulation  of  each 
peak  in  the  spectrum  individually  showed  that  the  low-frequency  peak  trans¬ 
formed  almost  identically  with  and  without  the  high-frequency  peak  present. 
The  energy  at  the  high-frequency  peak  decays  at  deeper  depths  when  the  low- 
frequency  peak  is  present.  Some  of  this  dissipation  takes  place  at  depths 
outside  what  is  traditionally  the  surf  zone. 

The  reason  for  the  interaction  between  the  two  peaks  is  not  understood. 

One  possibility  is  that  resonant  interaction  among  waves  transfers  energy  out 
of  the  higher-frequency  peak  to  even  higher  frequencies,  where  it  is  rapidly 
dissipated.  Another  possibility  is  that  the  bottom  friction  for  the  two-peak 
case  may  be  substantially  different.  However,  it  is  not  clear  why  all  the  loss 
should  be  at  the  high  frequencies.  A  third  possibility  is  a  shoaling  analog  of 
the  mechanism  proposed  by  Banner  and  Phillips  (1974)  in  which  the  high- 
frequency  waves  see  the  underlying  low-frequency  waves  in  terms  of  a  large- 
scale  flow  that  enhances  breaking  of  the  shorter  waves  with  no  apparent  effect 
on  the  longer  waves.  Which,  if  any,  of  these  explanations  is  correct  remains 
to  be  seen. 


Wave  Heights 


Figures  16  and  17  show  the  transformation  of  Hm  and  H^,  respectively, 
as  a  function  of  the  mean  water  depth  for  the  cases  with  dual  peaks  of  2.5  sec 
and  1.25  sec.  The  upper  three  curves  on  each  plot  are  the  high-energy  cases 
(total  initial  wave  height  of  15.2  cm)  with  three  different  distributions  of 
energy  between  the  two  peaks  (as  labeled  on  the  plot).  The  lower  three 
curves  are  the  low-energy  cases  (total  initial  wave  height  of  9.2  cm)  with  the 
same  three  distributions  of  energy  between  the  two  peaks.  Hw  tends  to  in¬ 
crease  as  the  water  depth  decreased  due  to  shoaling,  then  decreases  quickly 
when  the  depth  is  approximately  twice  the  height.  tends  to  stay  approxi¬ 
mately  constant  until  a  depth  of  about  twice  the  wave  height,  then  de¬ 
creases  quickly.  The  most  striking  feature  in  these  plots  is  the  difference  in 
H]a  due  to  the  different  distributions  of  energy  between  the  high-  and  low- 
frequency  peaks  (Figure  16).  The  maximum  difference  in  HU}  is  23  percent 
between  the  case  with  two-thirds  of  the  energy  in  the  low-frequency  peak  and 
the  case  with  one-third  of  the  energy  in  the  low-frequency  peak.  The  maxi¬ 
mum  difference  occurs  in  the  region  of  the  maximum  wave  height.  For  the 
dual  peaks  of  2.5  sec  and  1.75  sec,  the  different  distributions  of  energy  had  a 
smaller  effect.  The  maximum  difference  in  Hm  is  8  percent.  H^,  is  less 
influenced  by  the  distribution  of  energy  between  the  peaks.  The  maximum 
difference  in  is  13  percent  for  dual  peaks  of  2.5/1.25  sec  and  6  percent 
for  dual  peaks  of  2.5/1.75  sec. 
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Thus,  the  distribution  of  energy  between  two  peaks  can  have  a  significant 
impact  on  the  nearshore  wave  height  transformation.  The  impact  is  greater  on 
Hm  than  H^,.  The  proximity  of  the  peaks  also  affects  transformations.  Rela¬ 
tively  close  peaks  are  less  affected  by  different  distributions  of  energy  than 
more  separated  peaks. 


Summary 

Laboratory  simulation  of  two  wave  trains  shoaling  and  breaking  on  a  plane 
beach  indicates  that  there  are  strong  interactions  between  the  wave  trains. 
Simple  approaches  to  estimating  the  transformation  of  multiple  wave  trains, 
such  as  treating  the  wave  trains  independently  (and  superimposing  the  results) 
or  neglecting  one  of  the  wave  trains  can  lead  to  significant  errors  in  estimating 
wave  height  and  period  in  the  surf  zone. 

The  significant  wave  height  defined  from  an  analysis  of  individual  waves 
differed  by  more  than  20  percent  from  case  to  case  depending  on  the  locations 
and  relative  sizes  of  the  two  peaks  and  overall  energy  level.  This  is  of  engi¬ 
neering  significance  and  indicates  that  information  on  the  spectral  content  of 
the  wave  field  is  needed  to  make  accurate  predictions.  The  spectrally  based 
significant  wave  height  varied  by  about  10  percent  depending  on  location  and 
relative  size  of  the  two  peaks,  which  may  not  be  of  engineering  significance. 

Analysis  of  the  wave  spectra  indicates  that  the  low-frequency  peak  becomes 
dominant,  especially  in  shallower  water.  This  results  from  the  preferential 
loss  of  energy  in  the  high-frequency  peak.  Consequently,  methods  that  ignore 
the  lower  frequency  peak  or  superimpose  results  may  underestimate  the  wave 
period  in  shallow  water.  Mean  direction  and  directional  spreading  had  little 
influence  on  the  laboratory  results  for  these  plane-beach  experiments. 
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6  Conclusions 


This  report  summarizes  the  nearshore  wave  breaking  and  decay  research 
performed  under  the  Nearshore  Waves  and  Currents  work  unit.  The  work  has 
emphasized  wave  decay  through  the  surf  zone,  which  is  required  to  estimate 
nearshore  hydrodynamics  (longshore  and  cross-shore  currents  and  wave  set¬ 
up),  flooding  and  storm  damage,  and  nearshore  sediment  transport.  Incipient 
breaker  indices,  surf  zone  wave  decay  expressions,  spectral  shapes  in  the  surf 
zone,  wave  breaking  on  reefs,  and  decay  of  multiple  wave  trains  are  included 
in  the  report. 

Areas  for  future  research  include: 

a.  Wind  effects  on  wave  breaking  and  decay.  Existing  studies  on  wind 
effects  have  been  qualitative  and  have  not  provided  quantitative  rela¬ 
tionships. 

b.  Roughness  effects  on  wave  decay  over  reefs.  Field  measurements  of 
wave  decay  across  rough  coral  reefs  at  a  variety  of  locations  are  re¬ 
quired  to  quantify  wave  decay  due  to  bottom  friction. 

c.  Current  influence  on  wave  breaking.  The  impact  of  strong  longshore 
currents,  undertow,  and  ebb/flood  currents  (e.g.,  at  an  inlet)  on  depth- 
limited  wave  breaking  has  not  been  evaluated. 

d.  Detailed  modeling  of  the  energy  dissipation  in  breaking  waves.  To 
extend  wave  decay  modeling  beyond  empirical  methods,  the  details  of 
wave  breaking,  including  generation  and  dissipation  of  turbulence, 
must  be  included. 

e.  Incorporation  of  inffagravity  waves  into  wave  breaking  and  decay 
models.  Surf  beat  and  edge  waves  slowly  vary  the  mean  water  level 
and  the  mean  flow,  but  the  effect  of  these  long  waves  has  not  been 
included  in  wave  decay  modeling. 

f  Extending  nonlinear  wave  models  (e.g.,  Freilich  and  Guza  1984)  to 
the  surf  zone  with  continued  research  on  the  interaction  on  multiple 
wave  trains  in  the  surf  zone. 


Chapter  6  Conclusions 
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Appendix  A 
Agat  Laboratory 
Wave  Measurements 


A  laboratory  study  of  wave  decay  across  the  reef  at  Agat,  Guam,  was 
conducted  in  a  46-cm-wide  flume.  Starting  at  the  shoreward  end,  the  flume 
consisted  of  3.3  m  (covered  with  wave  absorber)  of  1/30  slope,  12.9  m  of 
reef  flat  (0  slope),  10.2  m  of  1/30  slope,  1.5  m  of  1/5  slope,  and  21.3  m  of 
flat  bottom.  The  waves  were  measured  at  six  locations  across  the  profile  at 
distances  of  38.0  m  (Gauge  1),  13.78  m  (Gauge  2),  12.04  m  (Gauge  3), 

9.66  m  (Gauge  4),  8.48  m  (Gauge  5),  and  6.25  m  (Gauge  6)  from  the 
shoreward  end  of  the  flume.  Table  A1  gives  the  zero-moment  wave  height 
H^,  peak  wave  period  Tp,  maximum  wave  height  H^,  and  the  mean  water 
level  t]  (relative  to  the  still-water  level)  for  69  test  cases.  The  still-water 
depth  h  in  the  deepest  portion  of  the  flume  is  also  given  in  Table  A1 .  Plots  of 
the  laboratory  measurements  and  wave  model  predictions  for  three  of  the  test 
cases  are  given  in  Figures  7-9  in  the  main  text. 
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Appendix  B 
Notation 


a  Empirical  variable  in  breaker  height  index  (Equation  6) 

b  Empirical  variable  in  breaker  height  index  (Equation  7) 

C  Wave  speed 

Cf  Bottom  friction  coefficient 

Ct  Wave  group  speed 

CtJ  Wave  group  speed  associated  with  the  stable  wave  height 
d  Mean  water  depth 

dh  Mean  breaker  depth 

E  Wave  energy 

E  Cg  s  Stable  energy  flux  associated  with  the  stable  wave  height 

E/W  Spectral  energy  for  FRF  spectrum 

/  Frequency 

fp  Peak  wave  frequency 

g  Gravitational  acceleration 

h  Still-water  depth 

hk  Still-water  breaker  depth 

H  Wave  height 

Ht  Breaking  wave  height 
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J7Mr  Maximum  wave  height 
H"  Zero-moment  wave  height 

Breaking  zero-moment  wave  height 
H",  Root-mean-square  wave  height 

Breaking  root-mean-square  wave  height 
Stable  wave  height 
Hc  Deepwater  wave  height 

Hta  Significant  wave  height 

Hl/10  Average  of  the  highest  1/10  wave  heights 

k  Wave  number 

L  Wavelength 

L,  Deepwater  wavelength 

n  Empirical  wave  decay  exponent  (Equations  IS  and  16) 

T  Wave  period 

Tp  Peak  wave  period 

x  Cross-shore  coordinate 

a.  Equilibrium  range  constant  (=0.0029)  (Equation  21) 

0  Beach  slope  in  radians 

/W  Spectral  wave  energy  parameter  for  FRF  spectrum 
y  Spectral  peakedness  factor 

7»  Breaker  depth  index 

r  Empirical  coefficient  defining  stable  wave  height-to-depth  ratio 

6  Energy  dissipation 

~r}  Mean  water  level 

6  Wave  direction  relative  to  shore  normal 
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k  Empirical  decay  coefficient  'Equation  18) 
£  Surf  similarity  parameter 

r  Constant  3. 14 159 

p  Density  of  water 

a  Spectral  width  parameter 

«  Radian  wave  frequency 

0*  Breaker  height  index 
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